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We consider cooperative processes (quantum spin chains and random walks) in one-dimensional 
fluctuating random and aperiodic environments characterized by fluctuating exponents uj > 0. At 
the critical point the random and aperiodic systems scale essentially anisotropically in a similar 
fashion: length (L) and time (t) scales are related as L ~ (rat) 1 / 1 ". Also some critical expo- 
nents, characterizing the singularities of average quantities, are found to be universal functions of 
uj, whereas some others do depend on details of the distribution of the disorder. In the off-critical 
region there is an important difference between the two types of environments: in aperiodic systems 
there are no extra (Griffiths)-singularities. 



I. INTRODUCTION 



Quenched disorder often has a strong effect on the co- 
operative properties of stochastic processes and strongly 
correlated systems, especially in one space dimension. 
For example the diffusion process in a one-dimensional 
environment (Sinai's walk) jl| becomes extremely slow, 
the mean-square displacement behaves like [X 2 (t)] av ~ 
In 4 1, in contrast to the linear i-dependence in the homo- 
geneous case. 

Ultraslow dynamics has recently also been observed 
in one-dimensional random quantum spin system 
The origin of the slow relaxation in these strongly corre- 
lated systems is again the presence of quenched disorder, 
and in particular the vicinity of a quantum critical point. 
Generally the presence of quenched disorder has a more 
pronounced effect on quantum phase transitions, which 
occur at zero temperature and are driven by quantum 
fluctuations, than on the so called classical phase transi- 
tions, which are driven by thermal fluctuations. 

Up to now, these two observations, namely the ultra- 
slow diffusion in general one-dimensional disordered envi- 
ronments and the so called "quantum" activated dynam- 
ics in random quantum spin chains, seemed to be unre- 
lated. One aim of the present paper is to demonstrate a 
intimate connection between both. Another issue is the 
question in how far these effects are also present in ape- 
riodic, i.e. non-random but nevertheless inhomogeneous 
environments. As we will see, relevant aperiodic systems 
bear a lot of similarities with completely random systems. 
However, also crucial differences exist, in particular in the 



so called off-critical regime, i.e. the Griffiths-McCoy 
region for the quantum spin chains and the anomalous 
diffusion regime of the random walk. It is simply not 
existent in aperiodic systems. A short account of the lat- 
ter issue (random versus relevant aperiodicity) has been 
given elsewhere Q. 

Although the results we report are of rather general va- 
lidity we have, for concreteness and for the lack of space, 
to confine ourselves to particular examples. So the pro- 
totype of a random quantum spin system is the random 
transverse Ising model (random TIM), which has received 
considerable interest recently ]7|-p^,p|,pd],p|p^r4]| . In one 
dimension the model is defined by the Hamiltonian 



H 



i+i 



hl<T? 



(1.1) 



in terms of the erf , erf Pauli matrices at site I. Here the 
exchange couplings J; and the transverse fields hi are 
quenched random variables, taken independently from 
distributions 7r( J) and p(h), respectively. The quantum 
model in (1.1) is closely related to a two-dimensional clas- 
sical Ising model with randomly layered couplings, which 
was first introduced and partially solved by McCoy and 
Wu (l|Jl|] (see also |T||j|]). The quantum control pa- 
rameter of the model is defined by 



5 = 



[ln/i] av -[lnJ] 

av 

var [In h] + var[ln J] 



(1.2) 



For S < (5 > 0) the system is in the ordered (disor- 
dered) phase and at S = there is a phase transition in 
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the system. One surprising observation is that at the crit- 
ical point some physical quantities are not self-averaging, 
thus the typical and the average values of those are dif- 
ferent. For these observables the critical properties are 
determined by the so called rare events, (which occur 
with vanishing probability) dominating the mean values. 

The average magnetization at the surface, m s and in 
the bulk, mi, of the system vanishes as a power law close 
to the critical point as m s ~ S@ 3 and nif, ~ 5@, respec- 
tively, with exponents 



f3 s = 1 and j3 = 



3 - V5 



(1.3) 



The average spin-spin autocorrelation function G(l) = 
[(erf of +r )]av involves the average correlation length £, 
which diverges at the critical point as £ oc |<5| -!y with 
the correlation length exponent 

v = 2, (1.4) 

which differs from the exponent of the typical correlation 
function G typ (l) = exp([ln(of af +r )] av }, which is 



Vtyp = 1 



(1.5) 



The time-dependent correlations of the model are very 
special, they differ completely from those in the pure sys- 
tem. The autocorrelation function 

G,(t) = [<*f (0)of (*)>W (1-6) 
at the critical point decays on a logarithmic scale 

Gi(t) ~ (lni)-" (S = 0), (1.7) 

where the decay exponent rj satisfies the scaling relation 
r\ = (3/v and i] s — (3 s /v for site I in the bulk and on the 
surface, respectively. Leaving the critical point in any di- 
rection one enters the Griffiths-McCoy regions, in which 
the connected autocorrelation function has a power law 
decay: 



(1.8) 



where the dynamical exponent z(S) is a continuous func- 
tion of the parameter 6. Close to the critical point it is 
z{5) w 1/25 . As a consequence of the power-law de- 
cay of the autocorrelations in the Griffiths -McCoy phase 
the magnetization is a singular function of the uniform 
magnetic field H x as m S i ng oc 

The above results about the random TIM are indepen- 
dent of the actual form of the probability distribution. It 
is often argued that the perturbation caused by the dis- 
order, with respect to the pure system, is connected to 
the fluctuating energy per spin: 



1 L 

A(L)/L=-J2(Ji-lJU)~L u 



(1.9) 



Here L is the linear size of the sequence of couplings Ji 
under consideration and uj — CL> ran d = 1/2 is the fluctu- 
ating or wandering exponent. The Harris criterion pcj ] 
makes a statement on behalf of the relevance/irrelevance 
of the perturbation of the critical behavior of the pure 
system by the above disorder. One compares the strength 
of the thermal fluctu ation s at the critical point with the 
fluctuating energy in fll.9|) on a length scale L identical to 
the correlation length. This yields in a one-dimensional 
system $ = 1 + — 1) for the cross-over exponent 
PHp2[ and indeed, for the random sequence the pertur- 
bation is relevant, since with v = 1 and u> = 1/2 one gets 
$ = l/2>0. 

It is known that there are non-random deterministic 
sequences, generated through substitutional rules, which 
have unbounded fluctuations, so that the corresponding 
fluctuating exponent is u> > 0. Having this similarity 
between random and aperiodic sequences in mind, one 
might ask the question whether or not the fluctuating 
exponent uj is the only quantity that determines the crit- 
ical behavior of systems with unbounded fluctuations in 
the couplings. 

As an example we consider the Rudin-Shapiro (RS) 
sequence p3j, which is built on four letters A,B,C and 
D with the substitutional rule: 

A -> AB , B -> AC , C^DB , D^DC. 

(1.10) 

Thus starting with a letter A one proceeds as: A — > 
AB -> ABAC -> ABACABDB -> etc., and one may 
assign different couplings to the different letters. The 
fluctuating exponent of the RS-sequence, which was orig- 
inally introduced to mimic random fluctuations, cjrs = 
1/2, i.e. just the same as for the random sequence. 

We note that up to now the critical behavior of ape- 
riodic systems has been studied mainly for such distri- 
butions that are non-relevant in the sense of the Harris 
criterion (i.e. uj < 0). Especially for marginal sequences 
(i.e. u> = 0) non-universal critical behavior coupling de- 
pendent anisotropy exponents have been found in exact 
calculations for the transverse field Ising chain [p4|-p7[ . 

In the relevant case, u) > 0, there are only a few exact 
results, obtai ned for one specific representation described 
after eq(l.lC) starting with the letter A. For this case the 
magnetization of finite RS-chains of length L at the criti- 
cal point has been shown to behave at the two end points 



77? 



: (RS) ~ exp(— const • VTj) 



(1.11) 



and m s (RS) = const > 0. Moreover, the lowest excita- 
tion energy at the critical point is given by 



(1.12) 



£i(£) ~ exp(— const • VL 



2=1 



We also note on recent Monte-Carlo simulations on the 
critical behavior of the Q = 8-state Potts model with 
layered aperiodic modulations Pgfl. 
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In the present paper we are going to study the aver- 
age quantities also for aperiodic chains. This means that 
for instance for the bulk magnetization we average over 
all sites of the system. This averaging process, however, 
could be performed in another way. One generates (hy- 
pothctically) an infinite sequence through substitution, 
cut out segments of length L starting at all different po- 
sitions of the (infinite) sequence and then average over 
all segments. Then, considering the example of surface 
magnetization one could expect that the two generic situ- 
ations in and below eq( 1.1 1| ) will appear with given prob- 
abilities, which arc then connected to the average critical 
behavior of the system. 

The structure of the paper is the following: in sec. II we 
present the basic notations of aperiodic sequences used 
in the paper. The free fermionic description of the TIM 
is given in sec. III. Results about the TIM both at the 
critical point and in the off-critical region are presented 
in sec. IV and V, respectively. We conclude with a final 
section on a parallel analysis of the behavior of a random 
walk in random or aperiodic environments. 



II. APERIODIC SUBSTITUTIONAL SEQUENCES 

In this paper we consider sequences generated via sub- 
stitution on a finite alphabet such that, in the case of 
two letters A and B, one substitutes A — > 5(A) and 
B — > 5(B). F or th e RS-sequence the explicit relations 
are given in eq(1.10). The properties of the sequence are 
governed by the substitutional matrix: 



_ "a "a 

~ T7 S(A) r7 5(B) 



where the matrix element nf (j) is the number of the let- 
ters i occurring in the substitution 5(j). The matrix ele- 
ments of M" contain the same numbers in the sequence 
after n iterations. 

If U r denotes the right eigenvector of M with eigen- 
values f2 r , the asymptotic density of i is given by 



(2.1) 



Ui(i) 

E,Ui(j) 



(2.2) 



where Ui is the eigenvector corresponding to the lead- 
ing eigenvector Oi. The length of the sequence after n 
iterations is related to the leading eigenvalue through 



(2.3) 



We consider the TIM in flLl| ) with a constant trans- 
verse field hi = h in an aperiodic environment. We assign 
to each letter in the sequence a coupling constant (for ex- 
ample Ja and Js). The cumula ted deviation A(L) from 
the average coupling [J] av in eq(1.9) scales with L as 



A(L) cx |n 2 |" oc L u 



(2.4) 



where 0,2 is the next-to-leading eigenvalue of the sub- 
stitution matrix M a nd oj is the fl uctu ating exp onent al- 
ready defined in fll.9|). Combining (2.4) and (2. 3) one gets 



In Iflal 



(2.5) 



The criti cal point h c of the system is obtained from the 
relation (1.2) with 5 = as 



X! pQ ln J i = ln h c 



(2.6) 



In the following we introduce a family of rele- 
vant sequences defined on A:-letters (fc-general sequence) 
Ai, A2, . . . , Afe with the substitutional rules: 



Ai 
A, 
A, 

A fc 



AiA 2 

A,-iA, + i 
Ai+iAi-a 
A fc A fe _i 



1 

k/2 



< i < 

< i < 



k/2 
k 



(2.7) 



The substitutional matrix of the sequence is given by 



I 1 1 
1 1 

1 1 



M : 



V 



with the eigenvalues 



1 1 
1 1/ 



(2.8) 



2cos[ T (r- 1)1 r = l,2,. 
k 



(2.9) 



Thus the fluctuating exponent uj from (2.5) for fc-general 
sequences is given by 



ln[2 cos(?r/fc)] 
la2 



(2.10) 



and is positive for k > 4. The leading right eigenvector 
of M is just Ui (i) = const, thus the asymptotic density 
of letters is given by 



= l/fc 



(2.11) 



We note that the /c-general sequence for k — 2 is just 
the Thue-Morse sequence, while for k — 4 we recover the 



RS-sequence mentioned in the introduction in (1.1C). 

To complete this section on the introduction of aperi- 
odic sequences we have to discuss the averaging process. 
As described in the introduction the average for aperiodic 
sequences is made in such a way that from a (hypothet- 
ical) infinite sequence segments of length L are cut out 
starting at all different points in the infinite sequence, 
these we call the realizations of length L of the aperi- 
odic sequence. Then we average physical observables over 
these realizations. The number of different realizations 
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R(L) of a random (binary) sequence grows exponentially 
with the size of the system, whereas in the aperiodic case 
have a much slower, namely linear increase: 



R(L) = ah 



(2.12) 



with a being some constant. For the RS-chain the num- 
ber of different realizations is less than 16L, and one can 
obtain the exact average value by the following proce- 
dure. Consider the RS-sequence in ( 1.10 ) generated from 
a letter A, take the first 4L sequels of length L, start- 
ing at positions 1,2, ...,4£ and take also their reflec- 
tion symmetric counterparts. Then repeat the procedure 
starting with a letter D. The averaging then should be 
performed over these 16L realizations. Therefore it is 
possible in numerical computations to do an exact aver- 
age for the RS-sequence, taking into account all different 
realizations and to obtain numerically exact results for 
relatively large (L < 512) aperiodic chains. For other 
values of k in the fc-general case the number of different 
realizations is still linear in L, however, no similar simple 
rule to generate all of them in linear computing time is 
known to us. 



III. FREE FERMIONIC REPRESENTATION OF 
THE TIM 



We consider the TIM in eq(l.l) on a finite chain of 
length L with free or fixed boundary conditions, thus 
we set Jl — 0. The simplest way to calculate various 
physical quantities of the model is to transform it into a 
fermionic representation according to Lieb, Schultz and 
Mattis |3(]] . A detailed description of the method and its 
application to different physical quantities can be found 
in ||], hereafter referred to as paper I. Here we briefly 
recapitulate the mai n re sults. 

The Hamiltonian (1.1) is expressed in terms of fcrmion 



creation and annihilation operators and n q , respec- 
tively: 



9=1 



(3.1) 



where the fermionic energies e q are identical to the non- 
negative eigenvalues of the 2L x 2L tridiagonal matrix T 
with eigenvectors V q : 



/Oh! 

hi Ji 

Ji h 2 

h 2 



Jl- 



Jl-i 


h L 



h L 




V a 



*,(2) 



(3.2) 



/ 



The magnetization profiles 



mi 



\of\0) 



(3.3) 



for symmetry braking boundary conditions (b.c.) can be 
expressed as I x I determinants, whose elements are de- 
fined in terms of th e co mponents \P 9 (i) and & q (i) of the 



eigenvectors V q in ( |3.2| ). If we fix the spin at one end of 
the chain, which is equivalent to setting hi, = 0, the local 
magnetization is given by 



free 



Hi Gu G12 
i?2 G21 G22 

Hi Gn G12 



Gu- 
G21- 



G, 



11-1 



(3.4) 



where 



Hj - $1(7) 



(3.5) 



The surface magnetization for fixed-spin b.c. m s — mi is 
given by the very simple formula 



L-l 



=1 j=i 



J, 



-1/2 



(3.6) 



which is exact for any finite system. 

The magnetization profile can be calculated analo- 
gously (see paper I), when both surface spins are fixed 
(++ and H — b.c.s). Finally , for a finite open chain the 
ground state expectation value of af in (3^) is identi- 
cally zero, due to symmetry. Here, as shown in Jl0| the 
matrix-element profile 



^frcc 



(ikflo) 



(3.7) 



is of importance, which is formally given by eq(3.4) with 
h L ^0. 

Here we note on a simple estimate for the excitation 
energy £i(L) in a open chain of length L p6| : 



4 



L-l 



Sx(L) ~ m s m s h L JJ 



(3.8) 



critical point has a similar scaling relation as in random 
systems: 



provided linii^oo e±(L)L = 0. Here m s and m s denote 
the finite size surface mag neti zations at both ends of 
the chain, as defined in eq(3.6) (for rn s simply replac- 
ing hj/Jj by /i/,_j/Ji_j in this equation). 

Next we consider the dynamical correlation functions 
of the system as a function of the imaginary time r: 

G,(r) = <0|(7f (rVf(0)|0) . (3.9) 
For surface spins this can be expressed in the simple form 



Gi(r) 



9 



|$ 9 (1)| 2 expire,) 



(3.10) 



whereas in the bulk G/ (t) can be expressed as a Pfaffian 
that can be evaluated via a determinant of / x I antisym- 
metric matrix (see paper I for details). 



IV. CRITICAL PROPERTIES 



(4.2) 



AE(L) ~ exp(-const • L u 



This form actually follows from the formula in eq( p.8[ ) for 
the energy gap of a sample with local order at the two 
ends, i.e. m s = 0(1) and m s = 0(1), where m s and m s 
are the surface magnetizations at the left and right end, 
respectively, to be calculated with eq(|3.6|). Then we have 
£ i ~ n'=i X ~ exp(— lt T \n{J/h)), where l tI ~ L u mea- 
sures the size of the transverse fluctuations a nd f n( J//i) is 
an averaged reduced coupling, from which (4.2) follows. 



To check the validity of the above scaling relation (4.2) 
we have investigated the probability distribution of the 
energy gap Pl(AE) at the critical point of the RS-chain 
and compared it with the same quantity for the random 
chain. 



In what follows we take for the Hamiltonian (LI) ho- 
mogeneous fields hi = h and two- valued couplings, say 
J i = A and J; = 1/A. For the fc-general sequence in (2.7) 
Ji = A for letters Ai with i < (k + l)/2 and J/ = 1/A for 
letters with i > (k + l)/2. For an odd k we ta ke J\ = 1 
for i = (k + l)/2. Thus for the RS sequence in (1.1CM we 
have J a = Jb = A and Jc — Jd — 1/A. Then from (2.6) 



and (2.11) for the critical point follows: 



In h r = 



(4.1) 



independently of k. Similarly, for the random model we 
take homogeneous fields and binary distribution of the 
couplings: J; = A and J; = 1/A, with the same prob- 
ability. The critical point in this case is also given by 
eq(O). 



A. The distribution of low energy excitations 

The basic features of the random transverse Ising chain 
at criticality have been described in the introduction. 
Here we stress again that many unusual properties of the 
random TIM are connected to the extremely broad distri- 
bution of the different physical quantities at the critical 
point. Concerning the energy gap AE in a finite system 
of size L the appropriate scaling variable is lnA-E/vL. 
The reason for this is that the system is essentially 
anisotropic at criticality, which means that times and 
length scale are related in an exponential rather than 
algebraic manner, c.f. eq( p..l2| ). 

For relevant aperiodic systems (i.e. those with a pos- 
itive fluctuatin g expo nent lo) one expects from scaling 
considerations plLp7 28 22] that the energy gap at the 



_ 0.1 

LU 

< 



a 



0.01 



0.001 



LU 

< 

a 



0.1 



0.01 



0.001 



o 




In AE / L 



1/2 




L= 8 
L = 16 
L = 32 
L = 64 
L= 128 



1 



1/2 



FIG. 



-In AE / L 

1. a) (Top) Scaling plot of the integrated probability 

for 

4, c.f. eq,fl4.3D, b) 



density fi(ln AE) versus the scaling variable (In AE/ L 1 ^ 2 
the RS sequence (exact average) with A 



(Bottom) The same as in a) for the the random sequence. 
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According to eq( |4.2j ) the appropriate scaling variable 
is hxAE/L 1 ' 2 for both systems. Indeed, as seen in Fig. 
la, the integrate probability distribution function 



Q L (lnAE) = 



In AE 



dyP L (y)^fl(\nAE/VL) (4.3) 



with Pl (In AE) = Pl(AE)AE, has a good data collapse 
using this reduced variable. Considering the same quan- 
tity for the random chain one can observe the same scal- 
ing behavior, as shown in Fig. lb. Thus we can conclude, 
that both the random and the RS-chains have logarithmi- 
cally broad distribution of the energy gaps at the critical 
point, from this fact one expects similar consequences for 
the critical behavior in the two systems. 



B. Surface magnetization 

For the TIM the surface magnetization represents per- 
haps the most simple order-parameter of the system, 
which in the fix-free b.c. is explicitly given by the simple 
formula in (3.6). For one single RS-chain, starting with 
the letter A this expression has been exactly evaluated 
at the critical point with the following results j28|. 

For A > 1, when the average couplings at the surface 
are stronger than in the bulk, the surface is ordered at 
the critical point, so that the surface magnetization is 
finite: 



lim m s (L, A, h = 1) 

L — >oo 



A 2 



1 



A 4 - A 2 + 1 



A > 1 (4.4) 



and approaches unity in the limit A — > oo. On the other 
hand for A < 1 , when the average couplings at the surface 
are weaker than in the bulk the critical surface magneti- 
zation vanishes as 



m s (L, A, h = 1) ~ exp(— const • V~L) 



(4.5) 



Proceeding by studying the critical point magnetizations 
in other sequels of the RS-chain one can notice that the 
above two examples are generic: in a sample the surface 
magnetization is either finite m s = 0(1), or it vanishes 
in the stretched exponential form as in eq(^j^). There- 
fore the average is dominated by the sample with finite 
surface magnetization ("rare events"), which occur with 
probability P r are(£) ^ L~^. Thus the critical surface 
magnetization is not self-averaging, it is determined by 
rare events and its scaling dimension x s m defined by the 
asymptotic relation [m s (L, h)] av ~ L~ Xm is just x s m = 7. 

In the following we make extended use of this obser- 
vation and calculated x s m exactly. Here we adopt the 
random walk picture of paper I. First we assign to each 
sample with a given realization of couplings a walk, which 
starts at the origin and makes the Z-th step +1 (—1) 
for a coupling J; = A (J; = 1/A). Second, we take 
the limit A — * 00, in which only those samples have 



non-vanishing surface magnetization, where the corre- 
sponding walk never visits sites with negative coordi- 
nates. Thus the proportion of rare events is given by 
the survival probability of the walk P raro (L) = P SUTV (L). 
Thus the third point in the study is to calculate the sur- 
viving probability of the walker. For the random chain 
the corresponding random walk is characterized by a sur- 
viving probability of P SU rv ~ L^ 1 / 2 , thus one gets the 
exact result: 



a4(rand) = 1/2 . 



(4.6) 



For the RS-chain we have performed the exact analysis 
of the average surface magnetization, the result of which 
is presented in the Appendix. The leading L-dependcncc 
of it is given by 

[m s (L, A - ex., h = l)] av = I (± + iV i- 1 / 2 

+ 0{L- 3/4 ) (4.7) 
from which the value of the scaling dimension 



< n (RS) = 1/2 



(4.8) 



follows. Thus we conclude that the surface magnetiza- 
tion scaling dimension is the same for the random and 
the RS-chain. 

In what follows we show that for relevant aperiodic se- 
quences x s m is a simple function of the fluctuating expo- 
nent w. We consider the scaling behavior of the surviving 
probability P SUIV (L) of the corresponding aperiodic walk, 
performing a discrete scale transformation, which corre- 
sponds to a substitutional step of the sequence. Then the 
length of the walk scales as L — > Lfii, whereas the trans- 
verse fluctuations Z tr scale like lt r — > Ztrl^l- Then the 
number of surviving walks with L steps, N(L), scales 
as N(L) -> N(LQt) = \Q 2 \N(L), since the number of 
these walks is proportional to the size of the transverse 
fluctuations. Remembering that the total number of 
different sequels is R(L) = aL the survival probability 
Psuiv(L) — N(L)/R(L) satisfies the scaling relation 



Psurv(^l^) 



■Psmv{L) 



P surv (L) ex L- (1 -^ 



(4.9) 



from which the value of the surface magnetization scaling 
dimension can be read off as 



= 1 



(4.10) 



We have seen by an exact analytical treatment that this 
relation is satisfied for the random and the RS-chain. For 
the first few me mber s of the family of fc-general sequences 
we verified eq(4.10) numerically and obtained x s m {k = 
5) = 0.3059(5) and x s m {k = 6) = 0.2076(4), which is in 
good agreement with the corresponding prediction from 
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(pLIOl) x s m (k = 5) = 0.3058 and < s (fc = 6) = 0.2075, 
respectively. 

Now we follow the analysis of paper I and calculate the 
correlation length critical exponent v from the S = In h 
dependence of the surface magnetization. In the scal- 
ing limit L > 1, 1 6 | -C 1 the surface magnetization can 
be written as [m s (L,S)] av — [m s (L, 0)] a . v rn s (6L 1 / 1 '). Ex- 
panding the scaling function into a Taylor series rn s (z) — 
1 + bz + 0(z 2 ) one obtains for the 5 correction to the sur- 
face magnetization: 



with 



[m s (L,6)] 

av 

[m s (L,0)] 

a 



e = i/v 



5L k 



(4.11) 



(4.12) 



This exponent can also be determined exactly in the 
A — > oo limit from random walk arguments. As shown in 
paper I the surface magnetization of rare events is given 
by: 



i {L,S) = {l + n)- 1 / 2 -5 



(n + 1) 3 /- 



+ 0(<5 2 ), (4.13) 



where the corresponding surviving walk returns n-times 
to its starting point after l\, I2, ■ ■ ■ , l n steps. The largest 
contribution to the coefficient of the term proportional to 
5 in (4.13) comes from those surviving walks, that have 



a large n, i.e. which visit the starting point frequently. 
For general aperiodic sequences n grows with the length 
of the walk as 



L a 



(4.14) 



and we call d s the surface fractal dimension of the sur- 
viving walks. One can check for the fc-general sequences 
that after two substitutional steps, when L — ► 4L, the 
nu mber of return points scales as n — > 2n, thus from 
eq(1.14) we get 



1/2 



(4.15) 



independent of the value k > 4. 

Next we ar e going to perform the average of the linear 
term in ( 4.13| ). Here we note that among the R(L) differ- 
ent samples there are D(L) = 0(n) that deliver the same 
dominant contribution: each of those has Nr(L) = 0(n) 
return points of characteristic length lj ~ l c hax (L) = 



0{L). Thus the average of the linear term in ( 4.13 ) grows 
like D(L) ■ l chal (L) ■ N R (L)/(R(L) ■ n 3 / 2 ) ~ L^Hence, 
comparing with ( [4.1 1| ) one gets the exponent relation: 



ds 
2 



(4.16) 



The r andom case ^(rand) = 2 is formally contained in 
(4.16) with d s = 0, since a surviving random walk re- 
turns n = 0(l)-times to the starting point. For the 
family of fc-general sequences d s = 1/2 for all values of 



fc, thus the corresponding correlation leng th exponent is 
v(k) = 4/(5 — 4wfe), with lo^ given in ( 2.10 ). In particular 
for the RS-sequence we get 



v{RS) = 4/3 



(4.17) 



We have checked this relation numerically by evaluating 
exactly the average surface magnetization up to L = 2 21 . 
The two point fits for the exponent O in (4.11) comparing 
systems of size 2 2l ~ 1 and 2 2l+1 are given in table I. Using 
standard sequence extrapolation techniques |31| on the 
data in Table I we got an estimate = 0.2501(2) which 
is in excel lent agreement with the scaling result (4.12) 
with ( Pjl ). 



I 


e 


2 


.19172235 


3 


.15877113 


4 


.16048277 


5 


.17656409 


6 


.19500760 


7 


.21068687 


8 


.22247716 


9 


.23086745 


10 


.23670532 



TABL E I: Numerical estimates of the exponent in 
(4.11) comparing numerically exact results on finite sys- 
tems of sizes 2 2l ~ x and 2 2l+1 . 



For the other members of the fc-general sequence the 
numerical estimates are v(k = 5) = 1.82(5) and v(k — 
6) = 2.23(5), what should be compared with the scal- 
ing predictions v(k = 5) = 1.799 and v(k = 6) = 2.186, 



obtained from eq(4.16) 



C. Magnetization profiles 

In a geometrically constrained finite system at the crit- 
ical point the appropriate way to describe the position 
dependent physical quantities, such as magnetization or 
energy density, to use density profiles rather then bulk 
and surface observables |32|. For two-dimensional clas- 
sical and one-dimensional quantum systems with homo- 
geneous couplings conformal invariance provides a useful 
tool to study various geometries. Let us consider a crit- 
ical system confined between two parallel plates, which 
are at large, but finite distance L apart, where the local 
densities cf>(r) vary with the distance I from the plates 
as a smooth function of l/L. According to conformal 
invariance lfc>3f 



m))ab = 



L . I 
— sm7T— 

7T L 



G ab (l/L) 



(4.18) 



where is the scaling dimension of the operator </> and 
the scaling function G a b(x) depends on the universality 
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class of the model and on the type of the boundary con- 
ditions to the left a and to the right b. With symmetric 
b.c. the scaling function is constant G aa = A. 

In two dimensions conformal invariance can also be 
used to predict the critical off-diagonal matrix elements 
profiles (<f>\ (f)(1) |0), where (</>| denotes the lowest ex- 
cited state leading to a non-vanishing matrix element 
(see These off-diagonal profiles give information 

about the surface and bulk critical behavior via finite size 
scaling. With symmetric b.c. one obtains for the profile 



(4.19) 



which involves also the surface scaling dimension x s d 



x m (RS) = 0.160(5) 



(4.20) 





r 


GD" 


sm7r— 




L 




l'/L 



FIG. 2. Scaling plot of the magnetization profile t71l(1) 
with symmetric (fixed at both ends) boundary condition 
(4.18) for the RS-sequence. The inse t shows a scaling plot of 
the off-diagonal matrix element (4.19) with symmetric (free at 
both ends) boundary condition. Here it A = 4, the bulk mag- 
netization exponent is used as a fit parameter: x m = 0.160(1) 
gives the best data collapse. 



For the random TIM model several profiles have been 
calculated in |l^,|| and they all follow very well the con- 
formal predictions. This coincidence of the numerical and 
the conformal results is quite surprising, since the ran- 
dom TIM is not conformally invariant due to strongly 
anisotropic scaling at the critical point. 

Here we consider the same problem for aperiodic chains 
and calculate the magnetization profiles for the RS-chain. 
In Fig. H we present the scaling plot of the diagonal 
and off-diagonal profiles with symmetric boundary condi- 
tions. Here we take the exact value for th e su rface mag- 
netization scaling dimension x s m from eq(4.8), whereas 



the bulk magnetization exponent, x m: is used as a fit 
parameter in order to obtain a good data collapse. As 
seen in Fig. |2| both profiles can be fitted with the same 
exponent 



which turned out to be independent of the parameter A 
and is different from the random chain value predicted 
by Fisher § to be x m (rand) = (3 - \/5)/4 = 0.191 . . .. 
The profiles, however, in analogy to the random case, 
follow very well the conformal predictions, both for the 
diagonal and off-diagonal profiles (4.18) and (4.19). This 
is an unexpected result, if we take into account that the 
relevantly aperiodic Ising chains are not conformally in- 
variant, due to anisotropic scaling at the critical point. 



D. Dynamical correlations 

Here we consider (imaginary) time depe ndent correla- 
tions of the same spin, as defined in eq(1.6). One expects 
different types of asymptotic behavior of the surface spins 
and of the bulk spins. First we consider the bulk auto- 
correlation function 



G(r) = [(o- x l/2 (t)*1 /2 ) 



'L/2l 



(4.21) 



and present a scaling consideration, where we essentially 
follow the steps of reasoning in the random case ||,^| . 

The autocorrelation function, like to the (local) mag- 
netization, is not self-averaging at the critical point: its 
average value is determined by the rare events, which 
occur with a probability P r and P r vanishes in the ther- 
modynamic limit. In the random quantum systems the 
disorder is strictly correlated along the time axis, conse- 
quently in the rare events with a local order, i.e. with 
a finite magnetization also the autocorrelations are non- 
vanishing. Under a scaling transformation, when lengths 
are rescaled as V = l/b, with b > 1 the probability of the 
rare events transforms as P' r = b~ Xm , like to the local 
magnetization. As we said above the same is true for the 
autocorrelation function: 



G(lnr) = b-^GQnT/b") 



5 = 



(4.22) 



where we have made use of the relation between relevant 
time t r and length £ at the cr itic al point, which follows 
from the scaling relation in eq(4.2). Setting the rescaling 
factor to b = (lnr) 1 /" we obtain: 



G(t) ~ (lnr)- 1 '"/" 



(4.23) 



For surface spins in ( p~22j ) and ( p~23| ) one should use the 
surface magnetization scaling dimension, x s rn . 

We have numerically calculated the magnetization au- 
tocorrelation function both at the surface and in the bulk 
for the RS-sequence. The results are depicted in Fig. ||. 
As can be seen the finite lattice results collapse onto one 
single curve (manifesting the absence of finite size effects) 
and the critical temporal decay happens on a logarithmic 
scale. The corresponding decay exponents are given by 
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x m /uo (xfjjui), for bulk (surfac e) co rrelations in agree- 
ment with the scaling result in ([4.23 ). 



law decay, 





FIG. 3. a) (Top): Bulk spin-spin autocorrelation function 
of the RS-sequence G 1 £/ 2 ( T ) ~ K C7 L/2(*) cr L/2)]av in imaginary 
time for various system sizes (and A = 4). Note that we 
have chosen L to be odd, so that L/2 denotes the central 
spin. In this plot with [G™/ 2 (T)]~ 1 ^ 2xm on lfn ear scale ver- 
sus r on a logarithmic scale the infinite system size limit is 
expected to lay on a straight line as indicated, b) (Bottom): 
Same as a) for the surface spin-spin autocorrelation function 
GT(r) = [(of(T)0f)]av in imaginary time. 

Next we consider the scaling behavior of the energy- 
energy autocorrelation function Gf(r) — [{af (r)af )] av . 
We note that af represents one part of the local energy 
operator, the other part of which — af af +1 — is related 
to it through duality, see paper I. Therefore the above 
autocorrelation function has essentially the same scaling 
behavior as the full energy density. 

For the random chain, as was shown in I, the critical 
energy autocorrelation function has an asymptotic power 



G?(r) 



(4.24) 



with the critical exponents 7/e(rand) 2.2 in the bulk 
and (rand) ~ 2.5 on the surface. For the aperiodic 
RS-chain the asymptotic decay is also consistent with a 
power law decay, as can be seen in Fig. |J, both in the 
bulk and at the surface of the system. The correspond- 
ing exponents 



r/ e (RS)«2.4 ?7|(RS)w3.4 



(4.25) 



are, however, different from those of the random chain. 




^ — 

a 





Q SL=17 




« 6L=33 




A A L=65 ! 




H hL=129 ^ 




o.n%- 3A 


ft 1 





100 
T 



a) (Top): Bulk energy-energy autocorrelation 
'L/2( r ) = [( cr l,/2( r ) cr i,/2)]av in imaginary time for 



FIG. 4. 

function Gf 

various system sizes (and A = 4) in a log- log plot. The straight 
line has slope —2.4, which yields our estimate for the exponent 
?7 e . b) (Bottom): Same as a) for the surface energy-energy 
autocorrelation function G\(r) = [{af (i")fi)]av in imaginary 
time. The straight line has slope —3.4, which yields our esti- 
mate for the exponent 
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V. OFF-CRITICAL PROPERTIES 



In random quantum systems there are Griffiths- 
McCoy singularities on the paramagnetic side of the crit- 
ical point, which result in a power law decay of the au- 
tocorrelation function G(t) ~ r -1 / 2 ^, where the dy- 
namical exponent z{8) is a continuous function of the 
control parameter 5. In the random TIM there is also a 
Griffiths-McCoy phase on the ferromagnetic side of the 
critical point and the values of the dynamical exponent 
in the two regions are related via duality, see paper I. 



1 



LU 

< 

a 



0.1 



0.01 



L= 16 
L = 32 
L = 64 
L = 128 



1 



3 4 5 
- In AE 




10 15 
- In AE 

FIG. 5. a) (Top) The integrated probability density 
fi(ln A_E) for the RS sequence (exact average) with A = 4 
slightly above the critical point (h — 1.5). The distribu- 
tion is chopped off at lnA_B mi 1 n (/i = 1.5) w 7.3. b) (Bot- 
tom) The same as in a) for the random chain with binary 
disorder, A = 4. Note the larger x-range as compared 
with the aperiodic sequence and the absence of a cut-off 
for large enough system size. The asymptotic form of the 
distribution is Q(x) ~ exp(x/z(8)), where x = lnA_B and 
l/z(ho = 1.5) = 0.40 (see paper I for details). 

The long time behavior of the average autocorrelation 
function G(t) is determined by the Laplace transform of 
the gap (ei) distribution function 



G(r) 



poo 

/ de P(e) cxp(— re) 
Jo 



(5.1) 



Thus the scaling properties of the low energy excitation 
are also connected to the above defined dynamical expo- 
nent z($): 



■(L,S) cx L~*W 



(5.2) 



in a finite system of length L. 

For relevantly aperiodic chains the same type of sce- 
nario, i.e. the existence of Griffiths-McCoy singularities, 
have been speculated |21|. In the following we are going 
to clarify this issue and study the distribution function 
of the energy gap £\ in the disordered phase (5 > 0) 
of the RS-chain. In Fig. | we compare the integrated 
gap distribution functions for the random and RS-chains. 
While the dat a fo r the random chain follow the scaling 
prediction in |jj) with S = 0.51n/i « 0.20, z w 2.5, 
the probability distribution of the RS-chain is chopped 
off: there is a L-independent cut-off at AE m i n (5). Con- 
sequently there is a relevant time scale in the problem 
t r ~ A.E~ in (<5) and the autocorrelation function has an 
exponential decay. The susceptibility and other physi- 
cal quantities are analytic in the whole disordered phase, 
thus there is no Griffiths-McCoy region in the RS-chain 
and we expect a similar behavior for any other aperiodic 
quantum spin chain. 

We can estimate the minimum energy gap AE m i n (S) 
as follows: We start with the formula for the excitation 
energy in ( |3.8| ) and consider a realization with surface 
order m s = 0(1) and fn s — 0(1), which is generally con- 
nected to the presence of a very small energy gap. Thus 
in the paramagnetic phase < S <C 1 we have 



A£ min ((5) ~ ei 



L-l 

n 



Ji 



exp(A5L - BIT) , (5.3) 



where A > 0, B > 0. The first term in the exponen- 
tial describes the average trend with S, whereas the sec- 
ond represents the largest possible fluctuation in the cou- 
plings among all aperiodic sequences of length L. It is im- 
portant to note that for a random distribution this second 
term could be proportional to L in rare events, implying 
that in the random case there is no minimum energy gap. 
For the aperiodic chains, however, due to the number of 
different realizations t hat i ncreases only linearly with L 
the fluctuating e nerg y (1.9) grows slower than L. Conse- 
quently from eq(5.3) one can derive a length scale 



L 



(5.4) 



that characterizes the most singular sample and the cor- 
responding minimum energy gap is then given by 



AE min (S) - exp (-const • S^^ 1 ^ 



(5.5) 



This relation is indeed satisfied for the RS-chain, 
since according to our numerical results in this case 
In AE min (6) ~ 1/6, as can be seen in Fig. [| 
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0.01 



0.1 



h-h 



o 



FIG. 6. The off-critical minimum energy gap of the 
RS-sequence lnAi5~ in (/i) versus the distance from the crit- 
ical point (h — ho ~ 5 for h — > ho) in the RS sequence for 
different s ystem sizes. The straight line has slope —1, as pre- 
dicted by (5.5). 



-p{t)=K-m) 



(6.3) 



Here 

P(t) = (P ifi (t),P i>1 (t), 



,Pi,L(t), p i,L+i(t)f (6.4) 



and the transition matrix is (M)j 



for i =/= j 



and (M)i 



J i,3 ■ 



while the initial condition is 
All physical properties of the model can 



be expressed in terms of (left and right) eigenvectors and 
eigenvalues of M, very much in the same way as the 
physics of the TIM is contained in th e ei genvectors and 
eigenvalues of the tridiagonal matrix ( |3.2| ). 

Here we consider first one quantity that gained consid- 
erable interest recently in related models for anomalous 
diffusion fl35| , [36"| ] : The persistence probability P pr (L,t), 
which is the probability that a walker starting at site 
z = l docs not cross its starting point until time t. Due 
to the absorbing sites at i = and i = L + 1 its long time 
limit p pr {L) = limt^oo P pr (L, t) is simply given by p7[ 



VI. RANDOM WALKS IN RANDOM AND 
APERIODIC ENVIRONMENTS 

As it was shown in paper I and being utilized in sec- 
tion IV/B in this paper there is a close relation be- 
tween the random quantum Ising spin chains and the one- 
dimensional random walk. Especially the scaling prop- 
erties of the surface magnetization and that of the low 
energy excitations of the random TIM can be obtained 
from the surviving properties of a one-dimensional ran- 
dom walk. Here we go further and emphasize a relation 
between the random TIM and the Id random walk in a 
random environment. 

To be specific, we characterize the one-dimensional 
random walk with the nearest neighbor hopping by the 
transition probabilities = w(i — > j) for a random 
walker to jump from site i to site j with 



Wi,i±l 





for 
for 



.71 = 1 

i| >i 



(6.1) 



Here we are particularly interested in the general case, 
in which the transition probabilities are not necessarily 
symmetric, i.e. 



+i,t 



(6.2) 



Moreover, the random walker is confined to a finite num- 
ber of sites i = 1, . . . , L. At the two ends of this interval, 
i.e. at i = and i = L + 1, we put absorbing walls, which 
is simply modeled by setting Wo,i = wl+i,l — (i.e. the 
walker cannot jump back into the system once landed on 
or L + 1). The time evolution of the probability dis- 
tribution of the walk Pi,j(t), which is the probability for 
the walker to be at time t on site j once started at time 
on site i, is fully determined by the Master-equation 



Ppr(L) = lim Pi,L+i(t) 

t—>oo 



i=l 3 = 1 W 3-3 + 1 



(6.5) 



Thus, as is shown in J37[ , there is a one-to-one relation 
between the persistence probabilit y (p.q ) and the surface 
magnetization m s (L) of the TIM 
correspondences 




with the following 



Wi,i+1 

m,i-i 

Ppr(L) 



hi 

ml(L) 



(6.6) 



Consequently similar relations hold for the average quan- 
tities, when the transition probabilities (or equivalently 
the fields and the couplings) follow the same random or 
aperiodic modulation. 

In the random case the critical point of the TI M c orre- 
spon ds to the Sinai walk |Q, and from eq(3.6), (4.6) and 
(|6.5| ) we have 



(6.7) 



For relevantly aperiodic environments which ar charac- 
terized by a wandering exponent to > we get from ( 4.10| ) 



^(aperiodic) (i)]av K L _(l- W ) _ ( g g) 

In the non-critical situation there is an average drift of 
the walk, which can be defined through eq(1.2) as 



few = 



[In W- 



[ln w,. 



varflni 



var[ln"> 



(6.9) 



where («;<_ ) stands for transition probabilities to the 
right (left), i.e. ii^j+i (wij-i). For <5rw 7^ the average 
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correlations denned on persistent walks are characterized 
by a correlation length 



3RW 



(6.10) 



with the exponents for the average given in eq(1.4) and 
eq(4.16) for the random and aperiodic environments, re- 
spectively. 

The dynamical properties of the random walk are dom- 
inated by the largest, non- vanishing eigenvalue A m of the 
transition matrix M. It can be shown B7n that for Am 



a formula similar to eq(3.8) holds, where hi and Ji have 
to be r eplaced by their random walk counterparts given 
in (6.6) and for the term m s mj stands the persistence 
pr obab ility p pr (L) . This implies straightforwardly, using 
cq([h^) that in a Sinai walk in the random and relevantly 
aperiodic case the characteristic time scales like 



(rand) 
char 



(id) 



exp (const • L") 



(6.11) 



with uj = 1/2 in the random case. As a further con- 
sequence, autocorrelations or return probabilities decay 
logarithmically, 



Po(t) 



i 



ln"(i) 



(6.12) 



with rj = 2 in the random case and rj = 1/lo in the 
relevantly aperiodic case. This follows from the scaling 
result P (t) 2 ~ [X 2 (t)]~^, where [X 2 (t)] m is the aver- 
age mean-square displacement, as already mentioned in 
the Introduction. 

Finally the Griffiths-McCoy phase of the random TIM 
is equivalent to the anomalous diffusion region of the ran- 
dom walk with S rw 7^ 0, in which case autocorrelations 
decay anomalously slow with an exponent 7(^rw) < 1 



J%(t)~f 



-7 (few) 



(6.13) 



that depends continuously on the drift parameter S^m 
and corresponds to the inverse dynamical exponent z(5) 
of the Griffiths-McCoy phase of the random TIM. In the 
two limiting cases we have z{8) = 1/2S as |<5| — * 
and z{5) — 1 as \5\ —too. In the latter case we approach 
the ballistic situation, when all the steps of the walk are 
made in the same direction, thus the average displace- 
ment has a linear time dependence. We can thus extend 
our dictionary (|6.6[) by 



A m 
7(<5rw) 



l/z(S) 



(6.14) 



For the aperiodic case with drift we remember the re- 
sults of section V, in particular eq(|l|) and conclude that 
a random walk in a relevantly aperiodic environment for 
<5rw 7^ does not exhibit a region of anomalous diffusion. 



VII. SUMMARY 



To summarize we have studied the effect of random 
and aperiodic environments on cooperative processes in 
one space dimension. We have shown that at the critical 
point, both for the transverse-field Ising model and for 
the diffusion process, the two types of inhomogencitics 
have quite similar consequences, which is based on the 
same type of distribution of the low energy excitations 
(large time scales). We have obtained - presumably ex- 
act - scaling relations, which connect the values of the 
surface magnetization exponent and that of the correla- 
tion length exponent with the known characteristics of 
the (random and aperiodic) environments. Besides the 
similarities between the critical properties of random and 
aperiodic models we have also observed several quantita- 
tive differences. For example some critical exponents are 
turned out to be environment dependent and — most 
noticeably — the Griffiths phase is absent for aperiodic 
models. 
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APPENDIX: AVERAGE CRITICAL SURFACE 
MAGNETIZATION FOR THE RUDIN-SHAPIRO 
CHAIN 



We consider the RS-chain in ( 1 -10| ) generated fro m a 
letter A and calculate the surface magnetization (3.6) 
for chains of length L = 2 2l+1 , I = 1,2, . . .. The chain 
starts at the 1st, 2nd,. . .,./Vth position of the original RS- 
chain, such that N = L ■ 2 2k , k — 1,2,..., and average 
over these N realizations. The average critical surface 
magnetization in the limit A — * oo is given by: 
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where 



m s (l, fc) = IjjVTxCi, fc) + + + H{k)[ Sl {l) + S 2 (l)} 



1 



n=l 
1 



+ 1 + n 



The asymptotic behavior of lim^oo m s (l, k) is given in eq(4.7) 



^ V2 fe -! +1+71 

ATi(Z, fc) = 18 • 4 fc " 2 + 9 • 2 k - 2 + 5(2' -1 - 1)(4* -1 + 2 k ~ 1 ) 
AT 2 (Z, fc) = A k - 2 (lO ■ 2 l ~ x + 4(7 - 1)) + (^5 ■ 2'" 1 + 2(Z - 1) - 2) • 2 

iV„(J, fc) = 2 2k - 1 + (4 fe - 1 - 2 fc - 1 )(2'- 1 - 1) 
H (fc) = 2 • 4 fe ~ 2 + 2 k - 2 



fc-2 



(Al) 



(A2) 

(A3) 

(A4) 
(A5) 

(A6) 
(A7) 
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